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Abstract 

We propose a new class of quantum computing algorithms which generalize many 
standard ones. The goal of our algorithms is to estimate probability distributions. 
Such estimates are useful in, for example, applications of Decision Theory and Artifi- 
cial Intelligence, where inferences are made based on uncertain knowledge. The class 
of algorithms that we propose is based on a construction method that generalizes 
a Fredkin-Toffoli (F-T) construction method used in the field of classical reversible 
computing. F-T showed how, given any binary deterministic circuit, one can con- 
struct another binary deterministic circuit which does the same calculations in a 
reversible manner. We show how, given any classical stochastic network (classical 
Bayesian net), one can construct a quantum network (quantum Bayesian net). By 
running this quantum Bayesian net on a quantum computer, one can calculate any 
conditional probability that one would be interested in calculating for the original 
classical Bayesian net. Thus, we generalize the F-T construction method so that it 
can be applied to any classical stochastic circuit, not just binary deterministic ones. 
We also show that, in certain situations, our class of algorithms can be combined with 
Grover's algorithm to great advantage. 
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1 Introduction 



In this paper, we use the language of classical Bayesian (CB) and quantum Bayesian 
(QB) nets[T]. The reader is expected to possess a rudimentary command of this 
language. 

We begin this paper with a review of various standard quantum computing 
algorithms; namely, those due to Deutsch-Jozsa[2], Simon[3|, Bernstein-Vaziranijl], 
and Groverj3]. We discuss these standard algorithms both in terms of qubit circuits 
(the conventional approach) and QB nets. Then we propose a class of quantum 
computing algorithms which generalizes the standard ones. 

Most standard quantum computing algorithms are designed for calculating 
deterministic or almost deterministic probability distributions. (By a deterministic 
probability distribution we mean one whose range is restricted to either zero or unit 
probabilities.) In contrast, our algorithms can also estimate more general probability 
distributions. Such estimates are useful in, for example, applications of Decision 
Theory and Artificial Intelligence, where inferences are made based on uncertain 
knowledge. 

Since some of the standard algorithms are contained in the class of algorithms 
that we propose, some algorithms in our class have a time-complexity advantage over 
the best classical algorithms for performing the same task. Even those algorithms 
in our class that have no complexity advantage might still be useful for nanoscale 
quantum computing because they are reversible and thus dissipate less power. Power 
dissipation is best minimized in nanoscale devices since it can lead to serious perfor- 
mance degradation. 

The class of algorithms that we propose in this paper is based on a construction 
method that generalizes a Fredkin-Toffoli (F-T) construction method |Hj used in the 
field of classical reversible computing. F-T showed in Refs.jO] how, given any binary 
gate / (i.e., a function / : {0, l} r — > {0, 1} S , for some integers r, s), one can construct 
another binary gate / such that / is a deterministic reversible extension (DRE) of 
/. / can be used to perform the same calculations as /, but in a reversible manner. 
Binary gates / and / can be represented as binary deterministic circuits. In this 
paper, we show how, given any CB net Af c , one can construct a QB net J\f® which is 
a "q-embedding" (q=quantum) of M c . By running on a quantum computer, one 
can calculate any conditional probability that one would be interested in calculating 
for the CB net J\f ' . Our method for constructing a q-embedding for a CB net is a 
generalization of the F-T method for constructing a DRE of a binary deterministic 
circuit. Thus, we generalize their method so that it applies to any classical stochastic 
circuit, not just binary deterministic ones. 

A quantum compiler [7j jS] can "compile" a unitary matrix; i.e., it can express 
the matrix as a SEO (sequence of elementary operations) that a quantum computer 
can understand. To run a QB net on a quantum computer, we need to replace the QB 
net by an equivalent SEO[2|. This can be done with the help of a quantum compiler. 
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Thus, the class of algorithms that we propose promises to be fertile ground for the 
use of quantum compilers. 

In certain cases, the probabilities that we wish to find are too small to be 
measurable by running J\f® on a quantum computer. However, we will show that 
sometimes it is possible to define a new QB net, call it M Q , that magnifies and 
makes measurable the probabilities that were unmeasurable using alone. We 
will refer to J\f® as Grover's Microscope for J\f®, because is closely related to 
Grover's algorithm, and it magnifies the probabilities found with J\[® . 



2 Notation and Other Preliminaries 

In this section, we will introduce certain notation that is used throughout the paper. 

We will use the word "ditto" as follows. If we say "A (ditto, X) is smaller 
than B (ditto, Y)", we mean "A is smaller than B" and "X is smaller than Y". 

Let Bool = {0, 1}. For integers a and b such that a < b, let Z a ^ = {a, a + 
l,a + 2,...6}. 

For any statement S, we define the truth function 9(S) to equal 1 if S is true 
and if S is false. For example, 6{x > 0) represents the unit step function and 
S(x,y) = 6(x = y) the Kronecker delta function. 

© will denote addition mod 2. For any integer x, x%2 will mean the remainder 
from dividing x by 2. For example, 4%2 = and 5%2 = 1. (This same % notation 
is used in the C programming language.) When speaking of bits with states and 1, 
we will often use an overbar to represent the opposite state: = 1,1 = 0. Note that 
if x, k G Bool then 

£(-l)**=l + (_l)* = 2<S(s,0). (1) 

k 

If x, y G Bool n , we will use x-y — Y.a=o x aUa, where the addition is normal, not mod 
2. 

Given x G Z 0tOO , let x = X)£L £ a 2 a , where x a G Bool for all a. Then we will 
denote the binary representation of x by bin(x) = (xq, Xi, £2, • • •)• Thus, bin a (x) = x a . 

On the other hand, given x = (xq, x±, x^, ■ ■ .) G Bool 00 , let x = Y^=o x a^ a ■ 
Then we will denote the decimal representation of x by dec(x) = x. 

We will use the symbol £), to denote a sum of whatever is on the right hand 
side of this symbol over those indices with a dot underneath them. For example, 
J2- /(?) b, c) = J2a,b /( a 5 b, c). Furthermore, J^aii wm denote a sum over all indices. If 
we wish to exclude a particular index from the summation, we will indicate this by a 
slash followed by the name of the index. For example, in Yl a ii/a,b we wish to exclude 
summation over a and b. Suppose / maps set S into the complex numbers. We will 
often use — to represent v^^ttt- Thus, num is shorthand for the numerator 

of the fraction. 
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We will underline random variables. P(a = a) = Pa{a) will denote the proba- 
bility that the random variable a assumes value a. P(a = a) will often be abbreviated 
by P(a) when no confusion will arise. Sa will denote the set of values which the 
random variable a may assume, and N a will denote the number of elements in S a . 
pd(B\A) will stand for the set of probability distributions P(-\-) such that P(b\a) > 
and J2b'eB P(b'\a) = 1 for all a G A and b G B. 

This paper will also utilize certain notation and nomenclature associated with 
classical and quantum Bayesian nets. For example, we will use (x.)a to denote {xi : 
i E A}. See Ref.[T] for a review of such notation. 

Hi — ^ | | ^ is the one bit Hadamard matrix. H^ B = Hf NB (the n-fold 

tensor product of Hi) is the Nb bit Hadamard matrix. We will also use Hi = -^H\ 
and H Nb = Hf NB = -j^=Hf NB . Note that {H^y = (-l) bb ' for b,b' G Bool, and 

(H NB ) i0 = {-lf 9 for b,V G BooI Nb . 

Any 2x2 matrix M which acts on bit a will be denoted by M(a). (We like to 
use lower case Greek letters for bit labels.) In this notation, a controlled-not (cnot) 
gate with control bit k and target bit r can be expressed as (7 x (r) n ^\ See Ref.jS] for 
more details about this notation. 

Let k = (k , k\, . . . , kn b _i) label iV^ bits. Assume all distinct. We 

will often use N$ = 2 Nb , where iV^ stands for number of bits and Ns for number of 
states. If \<p) Ki = \4>{. K i)) is a ket for qubit Ki, define \(f))n = \4>(k)) = IlS) -1 \4>i. K i))- 
For example, if 

|0>. = (J) (2) 

for all i, then 



Nb-1 / 1 \ / 1 \ (I 



|0>a= II |0>i«=( o)®lo)®"'®lo )=[ 1 .0.0.---.0] T - ( 3 ) 



i=0 



Likewise, if fi(«i) is an operator acting on qubit Ki, define Q(k) = n^o" 1 ^( K «)- For 
example, H\(k) = YI^q 1 H 1 (hi i ) is an Nb bit Hadamard matrix. 

Next, we will introduce some notation related to Pauli matrices. The Pauli 
matrices are given by: 

a * = ( i I ) ' = ( i ~o ) ' az = ( I -°i ) • (4) 

If \+ z ) and |— 2 ) represent the eigenvectors of a z with eigenvalues +1 and —1, respec- 
tively, then we define 

|0> = 1+,} = ( I ) , (5) 
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and 



|1> = !-*> = ( J J • (6) 



We denote the "number operator" by n. Thus 



and 



( I I \/ I l ~ a * (7\ 

«= oi J = I— -><— -I = — 2 — > ( ? ) 



n = l-n=[ . . =|+,)(+ z | = ^— 



n o / 

Since n and <r z are diagonal, it is easy to see that 

(-1)" = v. • (9) 
It is also useful to introduce symbols for the projectors with respect to |0) and |1); 

P 2 = |0)(0|=n, (10) 

Pf = |l)(l|=n. (11) 

Most of the definitions and results stated so far for a z have counterparts for 
a x and a y . The counterpart results can be easily proven by applying a rotation that 
interchanges the coordinate axes. Let w G {x,y, z}. If \+ w ) and |— w ) represent the 
eigenvectors of a w with eigenvalues +1 and —1, respectively, then we define 

\0 W ) = \+ w ) , (12) 

and 

IU = K> • (13) 

Let 

n w = 1 - n w = \+ w ){+ w \ = 1 + ° w ■ (15) 
As when w — z, one has 

(-1)"- = a w . (16) 

Let 
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Po=K)(O w \=n w , (17) 

and 

P™ = \l w )(l w \=n w . (18) 

In understanding Grover's algorithm, it is helpful to be aware of some simple 
properties of reflections on a plane. Suppose is a normalized (0'0 = 1) complex 
vector. Define the projection and reflection operators for by 

= 00+ , Rf = l- 2IL, . (19) 

Note that IT| = IF^. FigQ] shows that if x' = R^x, then x' is the reflection of x with 
respect to the plane perpendicular to 0. For example, R^cj) = — 0. 

X 

-2<|>«|> + x) 



Figure 1: Reflection with respect to plane perpendicular to 0. 

Some simple properties of Rs follows. = and R^R^ = R\ = 1. 

Since reflections are unitary matrices, a product of reflections is also a unitary matrix. 
Note that 




(-l) n * = e inTl * = 1 + (e inn * - 1) (20a) 

= 1 + II (e i7r - 1) = 1 - 211^ (20b) 

= iV (20c) 

(Eq. (|20b|) follows from the Taylor expansion of e" I " n< *.) 

If ei,e2,...,e n is an orthonormal basis for a vector space, IF = e;ej, and 
Ri = 1 — 211 i, then the product of the Ri in any order is —1. Indeed, 



R 1 R 2 ...R n = (l-2n 1 )(l-2n 2 )...(l-2n„) (21a) 

= i-2(nx + n 2 + ...n n ) (21b) 

= -1 . (21c) 
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Another property of reflection operators which is useful for understanding 
Grover's algorithm is the following. Let 




eo = ( J ) , ei = 



Now suppose that e[ is obtained by rotating e\ clockwise by an angle 9/2: 



_ ( cos(^/2) sin(fl/2) \ / sm((9/2) \ , , 

1 \-wx(9/2) cos(#/2) ) 1 ^ cos (6/2) J ' 1 ' 

e[ ~ ei for small 0. It is easy to check that the double reflection —R e / R eo is equivalent 
to a rotation (also clockwise) by 6: 

-Re>Re =( C ° S6 a Sin6 ) ■ (24) 

1 y -S1116 1 cos 9 J K ' 

(That these two successive reflections equal a rotation was to be expected, since 
the reflections are orthogonal matrices and a product of orthogonal matrices is itself 
orthogonal.) 

Above, we have considered plane reflections R^ acting on a complex vector 
space, but our formulas still hold true when R^ acts on a real instead of a complex 
vector space. In the case of real vector spaces, the Hermitian conjugate symbol f 
is replaced by the matrix transpose symbol T, and unitary matrices are replaced by 
orthogonal matrices. 



3 Some Standard Algorithms 

Next we will discuss several standard algorithms that are considered among the best 
that the quantum computation field has to offer at the present time. Later, we will 
try to generalize these standard algorithms. 

3.1 Deutsch-Jozsa Algorithm 

In this section we will discuss the D-J (Deutsch-Jozsa) algorithmj2]. We will do this 
first in terms of qubit circuits (the conventional approach), and then in terms of QB 
nets. 

Let k = (ko, kx, . . . , Kj\r B -i) label N B "control" bits and let r label a single 
"target" bit. Assume that r and all the Ki are distinct. We will denote the state of 
these bits in the preferred basis (the eigenvectors of a z ) by \x)%\y) T , where x G BooI Nb 
and y G Bool. Given a function / : BooI Nb — > Bool, define the unitary operator Q 

by 

Q = ^(^^(^^(^^^^(^^(^^(^^(r) , (25) 
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K 



G X H, 



Figure 2: Qubit circuit for D-J's algorithm. 



where n(n) = (n(n ), n(ni),..., n(rejv s -i))- The operation cr^ n ^ K ^(r), because it 
depends on /, is often called an "oracle" and each use of it is called a "query" . The 
right hand side of Eq.([25|) may be represented by the circuit diagram shown in Fig|5| 
The D-J algorithm consists of applying Q to an initial state |0)g|0) r of bits k and r, 
and then measuring the final state of these bits in the preferred basis. 

Fig|2]and the right hand side of Eq.([25[) are two equivalent ways of representing 
a particular SEO. There are infinitely many SEOs that yield Q. Figgis just one of 
them. In fact, the original D-J paper|2] gave a different SEO for Q, one with two 
queries instead of one. 

For X e BooI Nb ,Y e Bool, let 

|^o) = \X)z\Y) T , (26) 

and 

= Wi-i) for z = l,2,... , (27) 

where 

n x = H x [$H x (T)<T x (T) , (28) 
n 2 = a^(r), (29) 

and 

n 3 = n\ . (30) 

Then it is easy to show using simple identities (such as (#].)&,&' = (— l) bb ' , = 1, 
T = 0, and (-If = (-iy b for b, b' G Bool) that 
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N*) = ^+f E (-ir^-^'-^W| X% |y/ )r . (33 ) 
x,y,x',Y' 

Applying to the right hand side of Eq. (fH3j) and using the identity Eq.(JTJ) 

finally yields: 

{X',Y'\n\X,Y) = 6(Y',Y)-^ (-l) x < x '- x ) +¥T fW (34) 

xeBool N B 

for all X',X G BooI Nb and Y',Y G Bool. Thus, if the initial states of k and r are 
X_ = and H = 0, then the probability of obtaining X[ = X' for the final state of k 
is 



P(X'\X = Y = 0) = ^2\(x',Y'\n\x = o,Y = o)p 

y 
1 



^(_!)^'+/W|2 ( 35) 



Let Tbah the set of "balanced" functions, be the set of all / : BooI Nb — > Bool 
such that / maps exactly half of its domain to zero and half to one. Let T CO n-, the set 
of "constant" functions, be the set of all / : BooI Nb — > Bool such that / maps all its 
domain to zero or all of it to one. From Eq.(|35j). if X' = and / G Thai U Fcon, then 

P(X<=0|X=Z=0) = {J . ( 36 ) 

Now consider the QB net defined by FigJ3]and Tabled 
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nodes 


states 


amplitudes 


comments 


X 


X G BooI Nb 


5(X,0) 




Y 


Y G Bool 


6(Y,0) 




X 


x G BooI Nb 


(-1) x - x /V2 Nb 




y 


y G Bool 


(~l) yY /V2 




c 


c = (c x , c y ), c x G BooI Nb , c y G Bool 


8(c x ,x)5(c y ,y © /(x)) 




x! 


x 1 G BooI Nb 


<J(a/,ci) 




i 


y' G Bool 


^(l/. c s/) 




r 


X' G BooI Nb 


{-l) x '-* /V2 N * 




r 


Y' G Bool 







Table 1 



For this net, the amplitude A(x.) of net story x. is the product of all the terms 
in the third column of Tabled If 2£ = an d Y_ = 0, then the probability of obtaining 
X' = X' is 



P(X'\X = Y = 0) 



(37) 



num 

where A(x.) on the right hand side is evaluated at X = Y = 0. Substituting the 
value of A(x.) into Eq.([3T|) immediately yields Eq . (|33j) . 

Note that one can calculate the probability distribution Eq.([35j) by means of 
a CB net instead of a QB net. One can do this with the CB net defined by the graph 
X' -> Y' with: 
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nodes 


states 


probabilities 


comments 


x' 


X' e BooI Nb 


PAX') 




Y! 


Y' G Bool 


P ri x'(Y'\X') 





Table 2 



where Px 1 and Py'\x' are calculated from 

p x',Y'(X',Y') = \ (X',Y'\n\X = 0,Y = 0)| 2 . (38) 

We will say that the CB net defined by the graph X! ~* Yl and Table 121 is "q- 
embedded" in the QB net defined by FigJ3] and Table [T] In subsequent sections, we 
will say much more about q-embedding of CB nets. 



3.2 Simon's Algorithm 

In this section we will discuss Simon's algorithmjH]. We will do this first in terms of 
qubit circuits (the conventional approach), and then in terms of QB nets. 




Figure 4: Qubit circuit for Simon's algorithm. 



Simon's algorithm uses Nb "control" bits, just like the D-J algorithm. How- 
ever, it uses Nb target bits whereas the D-J algorithm uses only one. Simon's 
algorithm deals with a vector- valued function / : BooI Nb — > BooI Nb , whereas D- 
J's algorithm deals with a scalar-valued function / : BooI Nb — > Bool. Let K = 
(ko, kn b -i) label Nb "control" bits and let f = (r , 7~i, . . . , tjv b -i) label Nb 

"target" bits. Assume all n and Ki are distinct. We will denote the state of these 
bits in the preferred basis (the eigenvectors of a z ) by where x G BooI Nb and 

y G BooI Nb . Given a function / = (/o, /i, • • • , /jv b -i) where fi : BooI Nb — > Bool, 
define the unitary operator Q by 
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~Nb-- 

n 



T; 



HAk) 



(39) 



The operator Q for Simon's algorithm is analogous to the Q defined by Eq. for 
the D-J algorithm. The right hand side of Eq. (j39)l may be represented by the circuit 
diagram of Fig0] Simon's algorithm consists of applying Q given by Eq. (j39|) to an 
initial state |0),?|0)^ of bits k and f, and then measuring the final state of these bits 
in the preferred basis. One performs this routine several times. The measurement 
outcomes allow one to determine the period of the function / if / is of a special 
periodic type that will be specified later. 

Using the same techniques that we used to evaluate the matrix elements of Q 
for the D-J algorithm, one finds 



(X',Y'\Q\X,Y) 



2 N E 



E (-ir 



x-(X'-X) 



5(Y',Y®f(x)) 



(40) 



x£Bool N B 

for all X', Y', X,Y £ BooI Nb . If the initial states of k and f are X = and Y = 0, 
then the probability of obtaining X' = X' for the final state of k is 



P(X'\X = Y = 0) = J2\( X 'i Y 'M X = 0,Y = 0)| 2 

Y> 

= 7^EiE(-ir x '^/(^))r- (4i) 

^ yi x 

Now suppose Ts is the set of those functions / : BooI Nb — > BooI Nb such that 
/ is 2 to 1 (i.e., / maps exactly two domain points into each image point) and has 
a "period" A. By a period A, we mean a non-zero element of BooI Nb such that 
f{x) = f(x © A) for all x G BooI Nb . For any / G Fs and any y G BooI Nb , there 
exist exactly two elements of BooI Nb , call them X\ and X2, such that X\ = Xi © A and 
f(xi) = f(x 2 ) = y. Call f^iy) one of these x values, and call /~ 1 (y) © A the other. 
(The p subscript stands for "principal part", in analogy with Complex Analysis.) If 
/ G and /(/) is the image of /, then 



M y 6(f-\Y%x) + 6(f-\Y>)®A,x), if Y'el(f) 

0{Y,j{x))-<^ Q otherwise ' 1 ] 

Substituting this expression for S(Y', f(x)) into Eq.fjUJ) and using Eq.((TJ) yields 

P(X'\X = Y = 0) = ^5(X' • A, 0) . (43) 

To calculate the period A of /, run the experiment v times, measuring X 1 each time. 
Let X'{i) represent the ith measurement outcome. Then, for sufficiently large u, one 
can find A by solving the equations X'(l) • A = 0, X'(2) ■ A = 0, ... , X'(v) ■ A = 0. 
Now consider the QB net defined by FigEland Table El 
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Figure 5: QB net for Simon's algorithm. 



nodes 


states 


amplitudes 


comments 


X 


x e BooI Nb 


6(X,0) 




Y 


Y e BooI Nb 


8(Y,0) 




X 


x e BooI Nb 


(-1) x - x /V2 Nb 




c 


c = (c x , c y );c x , c y e BooI Nb 


6{c x ,x)6{c y ,Y ® f(x)) 




x! 


x 1 e BooI Nb 


6(x',c x ) 




r 


X' e BooI Nb 


(-1) x '- x '/V2 Nb 




r 


Y' e BooI Nb 


S(Y',c y ) 




Tab 


Le 3 



For this net, the amplitude A(x.) of net story x. is the product of all the terms 
in the third column of Table El If X — and Y_ = 0, then the probability of obtaining 
X' = X' is 



P(X'\X = Y = 0) 



Er Y,all/X'Y',X,Y A(x.)\ X =Y=0 



(44) 



J2x> num 

where A(x.) on the right hand side is evaluated at X = Y = 0. Substituting the 
value of A(x.) into Eq . (J44J) immediately yields Eq. (j4Tl) . 

It is possible to calculate the probability distribution Eq.(|41|) by means of a 
CB net instead of a QB net. One can do this with the CB net defined by the graph 
X' -> Y' with: 
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nodes 


states 


probabilities 


comments 


x' 


X' e BooI Nb 


PAX') 




Y! 


Y' G BooI Nb 


P ri x'(Y'\X') 





Table 4 



where Px> and Py'\x' are calculated from 

Px^Y'iX'X) = \(X',Y'\n\X = 0,Y = 0)| 2 . (45) 

We will say that the CB net defined by the graph X! —* Yl and Table 0] is "q- 
embedded" in the QB net defined by Fig|5] and Table El 



3.3 Bernstein- Vazirani Algorithm 

In this section we will discuss the B-V (Bernstein- Vazirani) algorithm^. 

To understand the B-V algorithm, it is helpful to first establish the following 
simple single qubit identities. First note that the single qubit Hadamard matrix 
rotates the Z-direction number operator into the X-direction number operator: 



Thus, 



[("I 



\n x -\b 



1 -1 
-1 1 

\bn x 



\bn z 



Hi 



Next note that a x exchanges the components of any vector it acts on: 



a 






a 



for any complex numbers a, 0. In particular, if b e Bool, then 

^10) = 1 6) • 



(46) 
(47) 

(48) 

(49) 



Now we are ready to discuss the B-V algorithm. Let k = (kq, k 1; . . . , «jv B -i) 
label Nb "control" bits and let r label a single "target" bit. Assume that r and all 
the Ki are distinct. We will denote the state of these bits in the preferred basis (the 
eigenvectors of a z ) by \x)z\y) T , where x G BooI Nb and y G Bool. For b G BooI Nb , 
define the unitary operator 



tot 



N B -1 

n a ^ n i 

i=0 



(50) 



The B-V algorithm is simply the following multi-qubit generalization of Eq. 
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wj|0)a = \b), . (51) 

That's all there is to B-V! 

Eq.([51|) can be represented by a qubit circuit consisting of a single wire for 
k, with a single node representing uo^. Eq.()51|) can also be represented by a QB net 
defined by the graph X_ — > X', with 



nodes 


states 


amplitudes 


comments 


X 


X=(X ,X l7 X NB ^)eBool N * 


6(X,0) 




r 


X' = (X' ,X' 1 ,X' NB _ 1 )eBool N » 


sHx'^) 





Table 5 



We should mention that it is common in the literature to dress up and obfus- 
cate Eq.fjHOJ) as follows. By virtue of Eq. fJTTj) . one can re-express cug as 

Nb ~ 1 N -1 

"i = II ^(^) 6l = ^i(^)(- 1 ) ElJ3 ° 6in * M #i(*0- (52) 

Some workers ascend to an even higher peak of obfuscation by adding a totally un- 
necessary target qubit. They define an operator, call it obtained by replacing the 
(— 1) in Eq.(JS2} by the operator a x {r) acting on a target qubit r: 

% = ^(^[^(r)]!^" 1 ^^^!^ • (53) 

At the beginning of the experiment, they put the target qubit in a state which is an 
eigenvector of u x {r) with eigenvalue —1. Thus, the obfuscated version of the B-V 
algorithm with a target qubit can be summarized by 

%K>r|0>* = wj|-x)r|0)a = \- x ) T \b) R . (54) 

We emphasize that for the B-V algorithm, the target qubit is a totally unnecessary 
affectation. 

So far we have given an unconventional presentation of the B-V algorithm. 
For completeness, we now give a conventional one. Define 

l^o) = |0)a|-*>T , (55) 

and 

= fiil^i-i) for i = 1,2,... , (56) 

where 

«i = #i(«) , (57) 



15 



and 



It follows that 



and 



n J = [ff,(T)]£S" ll ' n 'W 1 (58) 



n 3 = Qi . (59) 



l^l> = ^= E |5>«|-.>r, (60) 
l^> = ^E(-l) S *l^|-.>r, (61) 



i*> = ^u-ir^u-my)*\- x )r (62a) 



_i_^ ( _l)(?-^ g |_ x ) r (62b ) 
\b) R \- x ) T ■ (62c) 



To go from step (b) to step (c) of Eq.(JB2J), we used the orthogonality property given 
by Eq.flT). 

3.4 Grover's Algorithm 

In this section we will discuss Grover's algorithm [5] . 

Let k = (kq, ki, . . . , kat s _i) label A# bits. Assume all K,<i £1X6 distinct. We 
begin by defining the following N$- dimensional column vectors: 

\fi)x = V = »n s = -^=[1, 1, 1, 1 • • • , 1] T = ~^H Nb [1, 0, 0, 0, ... , 0] T , (63) 

|0)* = = [0,... ,0,0,1,0,0,. ..,0] T . (64) 

All components of are zero except for one predetermined component, located at 
position jtarg £ ^o,jv s -i, which equals one. We will refer to j tar9 as the target state 
(not to be confused with a target qubit). Note that we chose a special basis (or, 
equivalently, a special matrix representation) from the start. Note that = 
so /i and <p are nearly orthogonal for large Ng. It is also convenient to define the 
component-wise negation of 0: 
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|0™% = 0-' = [1,..., 1,1, 0,1,1,..., if. (65) 

Note that <p not is not normalized. 

Define projection and reflection operators for [i and </>: 

n M = |//)</x|, ^ = i - 2n„ = (-i) n - , (66) 

and 

n = |0)<0|, ^ = i-2n (/) = (-i) n ^. (67) 

Grover's algorithm can be summarized by the following equation|lUj 

(-R^Yfi « , (68) 

for some integer r to be determined, where "~" means approximation at large N$. 
Thus, starting with an Nb qubit system in a state /z, one applies the operator 
(—R^R^) consecutively r times, so that the N B qubit system ends in a state as 
close to <j> as possible. Measuring state in the special basis yields the target state 

Jtarg- 

Eq. can be represented by a qubit circuit consisting of a single wire for k, 
with r nodes, each representing —R^R^. Eq.fJHHJ) can also be represented by a QB 
net defined by a Markov chain graph X_o — > X_ l — > X 2 X r _ 1 , with 



nodes 


states 


amplitudes 


comments 




X G BooP s 


5(X ,0) 




X_i for i G Zi, r -i 


X G Boo^ 


Rfj.R<f>) Xj_i) 





Table 6 



To find the optimum number r of iterations, one can proceed as follows. 

First, notice that Eq. (j68|) describes a process which is entirely confined to 
the vector subspace spanned by /i and 0. Since /i and are not orthogonal, it is 
convenient to define an orthonormal basis eo,ei for the space span(fi,4>). Let 



knot 



^ = <P,e l = 7 ==. (69) 



Then 



fi=-=(e + ^N s -le 1 ). (70) 

FigEl portrays various vectors that arise in explaining Grover's algorithm. 

Since we plan to stay within the two dimensional vector space with orthonor- 
mal basis eo, e%, it is convenient to switch matrix representations. Within span(eo, ei), 
eo, ei can be represented more simply by: 
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Figure 6: Various vectors relevant to Grover's Algorithm. 



60 = ( I J ' 61 = ( 1 



If eo, ei are represented in this way, then 



1 



and 



— -R„iL 



iV 5 I -2v^V7^T iV 5 - 



Thus, 



R^Rff 



cos # sin 9 
— sin 6> cos 6 



where 



sin t7 = 



iV 5 VNs 

Eq . (|74j) is just Eq. (|24j) with e' x = /i and eo = 0. It follows that 



-Rfj,R<t>Y 



cos{rt>) smyrb 
— sin(r#) cos(r6 



and 



(—RfiR^y /J, 




sm(ro 
cos(rt 



v 7 ^ 
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We want the final state of the system to be parallel or anti-parallel to eo 
therefore, we want 



sin(r#) \ f ±1 
cos(r#) / 



(7? 



This will occur if 



r0«|(l + 2Jfc) , r « ^(1 + 2k)^N^ (79) 

for some integer k. 

Note that, in Grover's algorithm, the number of "queries" (calls to a unitary 
matrix that depends on 0) is far from unique. To illustrate this, let Q be a permuta- 
tion matrix that satisfies 

Q<t> = |0) = [1,0,0,. ..,0] T . (80) 
Since all the components of fi are equal, Q\i = \i. Thus 

(-R^Yfi = Q T (-R^R\Q)) r Qn = Q T (—R fJi R\ )) r fi . (81) 

Hence, it is possible to accomplish the full Grover transformation of /i with only a 
single query Q T . 

rotation by n/2. Let 



Grov 



U9 

-eieo + e e[ 



LnotiT , if inot\T 



12) 



Note that 



= 1 [— raot + (Ns — 1)0] ■ (83) 

^jVs(jV s - 1) 

From the point of view of quantum compiling, what Grover found is that the 7r/2 
rotation Ucrov is (approximately) equal to the r-fold product of —R^R^, where —R^R^ 
can be shown to have a SEO of low (polynomial in Nb) complexity. 
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Grover's algorithm has been modified in various, minor ways since it was first 
published. For example, Brassard et al. pointed out in Ref.^2] that the vector \l 
need not be the vector whose components are all equal. Other vectors \i will do just 
as well. Another modification of Grover's algorithm due to Younes-Miller |T3] adds 
an extra qubit to the original Nb qubits. Next we will discuss the Younes-Miller 
modification of Grover's algorithm, because it resembles a modification of Grover's 
algorithm that we will use in a future section. 

Let k = (k , K4,..., k Nb _i) label N B bits. Let r label a single bit. Assume 
r and all the Kj, are distinct. Let \i and denote the same N$ dimensional column 
vectors that we used in discussing the original Grover algorithm. In addition, define 
the following 2N$ dimensional column vectors: 



|# = |+.>r|A*>* = n ®^s= n > ( 84 ) 




10) = \~ X )r\4>)n = 7f ( -1 J ® ^ = 7f ( -0 J • 

Note that (<t>\p) = ^2N S ' so ^ an< ^ & are near ^y orthogonal for large Ns- Define 
projection and reflection operators for in the usual way: 

11^ = |0) (0| , R} = 1 - 21^ = 1 - 2IL w M hx)T . (86) 
R^ can be re-expressed as 

R 4, = 1 + U \<f>)A a x( T ) - 1) = ex P[ n |*> K lna x (r)] = 

= W^- (87) 

Define projection and reflection operators for jl in the usual way: 

II A = \p){p,\ , R fi = 1 - 2II A = 1 - 2U Mx U 1+z)t . (88) 
Rjx can be re-expressed as 

R fi = H^k) (l - 2n| 0)B n|o) T ) H^k) . (89) 

In analogy with the original Grover's algorithm, the Younes-Miller version can be 
summarized by 

(-RpRiYlii) « |0) , (90) 

for some integer r to be determined, where means approximation at large iVg. 
Thus, starting with an Nb + 1 qubit system in a state /2, one applies the operator 
(—R^R^) consecutively r times, so that the final state of the Nb + 1 qubit system 
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ends in a state as close to <ft as possible. Measuring state <f> in the special basis yields 
the target state jtarg- 

To find the optimum number r of iterations, one can proceed as follows. 

First, notice that Eq. (|90j) describes a process which is entirely confined to 
the vector subspace spanned by ft and 0. Since ft and <fi are not orthogonal, it is 
convenient to define an orthonormal basis eo, ei for the space span(ft,4>). Let 



and 



e = <f) 



1 



J_( 

y/2 [ -4 



ei = J?\fi- (/f e o)e ] , 
where K is chosen so that e\ = 1. It is easy to show that 

K — [ft — (ft ■ e )e 



\ N s 



Thus, 



1 



ei 



Furthermore, 



1 



e 



2iV. 



lei 



(91) 



(92) 



(93) 



(94) 



(95) 



Fig0 portrays various vectors that arise in explaining Younes' version of 
Grover's algorithm. 




Figure 7: Various vectors relevant to Younes' version of Grover's Algorithm. 



Since we plan to stay within the two dimensional vector space with orthonor- 
mal basis eo, ei, it is convenient to switch matrix representations. Within span(eo, ei), 
e , t\ can be represented more simply by: 
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eo = ( o J , ex = ^ I J . (96) 
If eo, ei are represented this way, then 




^m~s \ V2N S - 1 



and 



Thus, 




(97) 



_ 1 / Ns-1 #R \ (qR] 

- R * R i> - w s [ N S -i ■ (98) 



-R.R.-l cosd sin M fgg) 



where 



sin^ = ^ wJ — . (100) 

N s y N s V } 

A comparison of Eq. (|72jl (for the original Grover's algorithm) and Eq. ()97j) (for 
Younes's version of Grover's algorithm) reveals that for the purpose of finding the 
optimal number r of iterations, Younes' algorithm is the same as Grover's algorithm if 
one replaces N$ in Grover's algorithm by 2N$- This comes from the fact that Younes' 
algorithm uses iV^ + 1 bits whereas Grover's uses A^^. 



4 Generalization of Standard Algorithms, 
a list of Desiderata 

So far we have analyzed several standard quantum computing algorithms, namely 
those attributed to Deutsch-Jozsa, Bernstein- Vazirani, Simon and Grover. (Two 
other standard algorithm's that we didn't analyze are Shor's algorithm [Tij and the 
algorithm for Teleportation J3].) In this section, we will try to point out those features 
of the standard algorithms that would be, in our opinion, fruitful to generalize. Bear 
in mind that generalizations are seldom unique, but some are more natural, fruitful 
and far-reaching than others. 



(a) Allow more complicated graph topologies 



The standard algorithms discussed here can all be represented by QB nets 
with trivial topologies such as 2 body scattering graphs or Markov chains. However, 
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other important quantum algorithms, such as the one for Teleportation^l, can be 
represented by QB nets with more complicated graph topologies (e.g., with loops). 

(b) Estimate more general probability distributions 

The goal of most standard algorithms is to estimate a deterministic probability 
distribution. However, estimating non-deterministic ones is also very useful. Such 
estimates are useful in, for example, applications of Decision Theory and Artificial 
Intelligence, where inferences are made based on uncertain knowledge. 

(c) Allow multiple runs and the rejection of some 

If one is estimating a non-deterministic probability distribution, it will be 
necessary to do multiple runs. It may also be necessary to allow rejection of runs. 
Obviously, the number of rejected runs is best kept as small as possible. 

(d) Allow more general measurements 

Suppose x is a node of a QB net. Let S x be the set of its possible states. 
We will say that node x has been measured if during the experiment which the QB 
net describes, a measurement is performed that restricts the possible states of x to a 
proper subset S' x of S x . When x is an internal (ditto, external) node of the QB net, 
we will refer to its measurement as an internal (ditto, external) measurement. 

The standard algorithms discussed here use external but no internal measure- 
ments. However, other important quantum algorithms, such as the one for Telepor- 
tation, do use internal ones. 

5 Q-Embeddings 

The remainder of this paper will be devoted to discussing a class of algorithms which 
generalizes some standard algorithms and achieves some of the desiderata given in 
the previous section. Our algorithms are based on the idea that, given a CB net, 
one can always embed it in a QB net. Simple examples of such q-embeddings have 
already been given in the sections dealing with standard algorithms. 

We start by defining some terminology that will be useful. 

A probability matrix P(y\x) is a rectangular (not necessarily square) matrix 
with row index y G S y and column index x G S x such that P(y\x) > for all 
x, y, and Y^yP{y\ x ) — 1 f° r & H x - The set of all probability matrices P(y\x) where 
x G Sx and y G S y will be denoted by pd(S y \S x ) (pd = probability distribution). 
A probability matrix is assigned to each node of a CB net. A probability matrix 
P(y\x) is deterministic if for each column x, there exists a single row y, call it 
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y(x), such that P{y\x) = 5(y(x),y). Any map / : S^—* S y uniquely specifies (and is 
uniquely specified) by the deterministic probability matrix P with matrix elements 
P(y\x) = 5(y, f(x)) for all x G Sx and y G S y . We will often talk about a map / and 
its associated probability matrix P(y\x) as if they were the same thing. 

Given two matrices A and B of the same dimensions, their Hadamard prod- 
uct C = AQB is defined by C id = A^B^ for all We will call HAS (A) = AQ A ] 
the Hadamard Absolute Square (HAS) of matrix A. If U is a unitary matrix, 
then HAS (£7) is a probability matrix. For example, for any angle 9, 



HAS( 



Another example is 



cos 9 sin 9 
- sin 9 cos 9 



HAS(#i 



cos 2 9 sin 2 9 
sin 2 9 cos 2 9 



1 1 
1 1 



(101) 



(102) 



A CB net M c is the HAS of QB net M Q if M Q and M c have the same 
graph, and their node matrices are related as follows. For each node x i7 if 74[xj|(x.)rJ 
is the amplitude of node x { in J\f Q , and P[xj|(x.)rJ is the probability of node x { in 
M c ', then |A[xj|(x.) ri ]| 2 = P[xj|(x.) r J. In such a case, we will write HAS(A/" Q ) = M c . 

A unitary matrix A(y, x\x, y) (with rows labelled by y, x and columns by x, y) 
is a q-embedding of probability matrix P{y\x) if 



E \Mv, 



x\x, y 



P{y\x) 



(103) 



for all possible values of y and x. (the "q" in "q-embedding" stands for "quantum"). 
We say y is a source index and x is a sink index. We also refer to x and y collectively 
as ancilla indices. Note that any unitary matrix is a q-embedding of its HAS. Indeed, 
in this case Eq. ()103J) is satisfied with the indices x and y each ranging over a single 
value (i.e., x and y are fixed). If a q-embedding satisfies A(y,x\x, y) G Bool for 
all y, x, x, y, we say that it is a deterministic q-embedding or a deterministic 
reversible extension (DRE) of its probability matrix (note that its probability 
matrix must also be deterministic). By an extension of a matrix we mean adding 
extra rows and/or columns to it. General q-embeddings use the square root of the 
entries of the original probability matrix so they are not simply extensions of the 
original matrix; they are, however, reversible since they are unitary matrices. 
Given a QB net Af Q , let 



P[(x.) 



E 

O0r o - 



A{x.) 



(104) 



On the right hand side of Eq. (jl04j) . A(x.) is the amplitude of story (x.), Tq is the set 
of indices of all the nodes of Af Q , and L is the set of indices of all leaf (aka external) 
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nodes of Af Q . We say Af Q is a q-embedding of CB net Af c if P[(x.)l] denned by 
Eq.flnHJ) satisfies 

P[(x.) rc } = Y,P[(x.) L ], (105) 

Li 

where Li C and is the set of indices of all nodes of Af c ' . Thus, the probability 
distribution associated with all nodes oiM c can be obtained from the probability dis- 
tribution associated with the external nodes of Af®. Some examples of q-embeddings 
of CB nets have already been given during our discussion of standard algorithms. 
More examples will be given in subsequent sections. 

For some positive integers r and s, we will say a map / : BooV — > Bool s 
is a binary gate from r to s bits. / uniquely specifies (and is uniquely specified) 
by the deterministic probability matrix with entries P(y\x) = S(f(x),y), where x = 
(xo, xi, . . . , x r -i) G BooV and y = (yo, yi, ■ ■ ■ , y s -i) £ BooV. If / is an invertible 
map, we will say that the gate is reversible. For example, the AND gate which 
takes (xi,x ) — > y with y = x Xi is a binary gate. So are the OR and NOT gates. 
Out of these 3 gates, only the NOT gate is reversible. 

Another example of a reversible binary gate is the Toffoli gate[6]. It maps 3 
bits into 3 bits as follows: 

Vo = T (x) = x , 

y 1 = T 1 (x) = x 1 , (106) 
y 2 = T 2 (x) = x 2 © xqXx. 

The Toffoli gate can also be defined as the following deterministic probability matrix 

P(y\x) = 6(y,T(x)) = 5(y 2 ,x 2 © x x 1 )6(y 1 ,x 1 )6(y ,xo) . (107) 

Consider 3 bits labelled 0, 1, and 2, and suppose the ith bit changes value from x\ to 
yi. Then bits and 1 do not change whereas bit 2 flips iff the product xqX\ equals 
one. Thus, the probability matrix with entries given by Eq. ()107|) is simply a doubly 
controlled not: 

[P(y\x)]=o x (2) n ^ . (108) 

It is convenient to use the term Toffoli gate to refer not only to the gate defined by 
Eq. (jl07|) . but also to the 3 other gates that one obtains by replacing x Xi in Eq. (jl07|) 
by XoX~i, or XoXi, or XoX±. This corresponds to replacing n(l)n(0) in Eq. (jl()8j) by 
n(l)n(0), or n(l)n(0), or n(l)n(0). FiglHl shows the 4 doubly-controlled nots that we 
call Toffoli gates as well as the circuit diagrams usually used to represent them. 

5.1 Q-Embedding of Probability Matrices 

In this section we will first give some examples of q-embeddings of probability matri- 
ces. Then we will show that any probability matrix has a q-embedding. 
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n(l)n(0) 
2 1 



u u u 



n(l)n(0) 

^(2) 

2 1 



U u u 



-6 



n(l) n(0) 

a (2) 

2 1 



n(l) n(0) 
2 1 



Figure 8: Four different kinds of Toffoli gates. 0,1,2 are bit labels. 



Any unitary matrix is a q-embedding of its HAS, but such q-embeddings are 
trivial in the sense that they have no ancilla indices. 

As first shown in Refs.jH], the Toffoli gates can be used to build q-embeddings 
(in fact, DREs) of the elementary binary gates AND, XOR, NOT, FANOUT. See 
FigEl Let x = (xq, Xi,X2) € Bool 3 and y = (yo, y%, y 2 ) £ Bool 3 . For the AND gate, 



E |<!/k.(2) 



n(l)n(0) I 



X 2 = 0, Xi, Xq) = 8(y 2 , XiXq) 



yi.3/0 



For the FANOUT gate, 



E I (V l^-(2) n(1)n(0) |z 2 = 0, X! = 0, z ) = <% 2 , zoWj/o, x ) . 



yi 



For the XOR gate, 



E |d/k.(2) 



n(l)n(0)| 



x 2 ,xi,x = 0) = <%2,£2 ©Si) • 



yi,yo 



For the NOT gate, 



(109a) 



(109b) 



(109c) 



E |<i/M2) 



n(l)n(0) I 



x 2 ,xi = 0,x = 0) = <J(y 2) ^2 © 1) = S(y 2 ,xJ) . (109d) 



3/1,2/0 



Note that the NOT gate is just a x , which is a DRE of itself. Eq. (|109djl gives a 
different DRE of a x . In the left hand side of Eqs. (jl09}) . the Xi indices that are set to 
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AND 
X, X, 



y- 



X i X o 



x, X, 



FANOUT 
X, 



X 6 O 

T T " 



X n x r 



XOR 

X, X, 



-6 



x 2 ©x 1 X, 



NOT 
X, 



X 6 6 

I T ▼ 



x o 



Figure 9: How to express elementary gates in terms of Toffoli gates, xq, x±, X2 G Bool 
are bit values. 



zero are called source indices, and the indices that are summed over are called 
sink indices. Sink and source indices are collectively called ancilla indices. 

Next we will prove that any probability matrix has a q-embedding. Suppose 
that we are given a probability matrix P(y\x) where x G S x and y G S y . Let N x 
(ditto, N y ) denote the number of elements in Sx (ditto, S y ). Let for x G S^ be 
any orthonormal basis of the complex N% dimensional vector space. The components 
of £^ will be denoted by £g , where x G S%- If the ^ x ^s are the standard basis, then 
= 8(x, x). Define matrix A by 



A(y,x\x,y) = \ V P (^) & if ^ = . (110) 

( obtained by Gram-Schmidt method if y ^ 

To understand the last equation, consider FigllOl In that figure we have 
assumed for definiteness that S x = {0, 1, 2} and S y = {0, 1, 2, 3}. The shaded (ditto, 
unshaded) columns have y ^ (ditto, y = 0). It is easy to see that the unshaded 
columns are orthonormal because the vectors £ x are orthonormal and J2 y P{y\ x ) — 
1. Since the unshaded columns are orthonormal, one can use the Gram-Schmidt 
method [TU] to fill the shaded columns so that all the columns of A are orthonormal 
and therefore A is unitary. Note that by virtue of Eq. (|110j) . 

E \A{y,x\x,y = 0)| 2 £P(y|s)d"M" = P(v\x) (HI) 
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Figure 10: How to construct a q-embedding of any probability matrix. 



so that the A defined by Ea. ffTTUl) does indeed satisfy Ea. (fTU3jl .[Tfl] 

Note that the matrix A defined by Eq. fjllOj) has dimensions N%_N y x N^N y . It 
is sometimes possible to find a smaller q-embedding of an N y x Nx probability matrix 
P(y\x). For example, a x is a q-embedding of itself. As a less trivial example, suppose 

P(y\x u x 2 ) = 5{y,x 1 ®x 2 ) , (112) 
for y,x±,x 2 G Bool. Then define 

f_iyie 

A(y,e\x 1 ,x 2 ) = -j= 5(y,x 1 ®x 2 ), (113) 
for y,e,xi,x 2 G Bool. It is easy to check that matrix A is unitary. Furthermore, 

^2\A(y,e\x h x 2 )\ 2 = -J2S(y,x 1 ®x 2 ) = 5(y,x l ®x 2 ) . (114) 
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5.2 Q-Embedding of CB Nets 

As we've said before, F-T showed in Refs.jS] how, given any binary gate /, one 
can construct another binary gate / such that / is a DRE of /. Their method for 
constructing / is to first represent / as a binary deterministic circuit composed of 
elementary gates (AND, XOR, NOT, FANOUT), and then to modify the circuit by 
replacing each of its gates by a DRE of it. The desired gate / is then specified by the 
modified circuit. 

In this section we will show how, given any CB net Af c , one can construct a 
QB net J\f® which is a q-embedding of M . So far we've shown how to construct a 
q-embedding for any probability matrix. Now remember that each node of M c has 
a probability matrix assigned to it. The main step in constructing a q-embedding 
of M c is to replace each node matrix of M with a q-embedding of it. Thus, our 
method for constructing a q-embedding of a CB net is a generalization of the F-T 
method for constructing a DRE of a binary deterministic circuit. We generalize their 
method so that it can be applied to any classical stochastic circuit, not just binary 
deterministic ones. 

Before describing our construction method, we need some definitions. We say 
a node m is a marginalizer node if it has a single input arrow and a single output 
arrow. Furthermore, the parent node of m, call it x, has states x = (xi,x%, . . . ,x n ), 
where X{ G S^. for each % G Z\ >n . Furthermore, for some particular integer i G Zi jTl , 
the set of possible states of m is 5 m = S^, , and the node matrix of m is P(m = 
m\x = x) = S(m, x io ). 

Let M c be a CB net for which we want to obtain a q-embedding. Our con- 
struction has two steps: 

(Step 1) Add marginalizer nodes. 

More specifically, replace Af c by a modified CB net Af c mo d obtained as follows. 
For each node x of M c , add a marginalizer node between x and every child of x. If 
x has no children, add a child to it. 

As an example of this step, consider the net M c ("two body scattering net") 
defined by FigJTT]and Table 



nodes 


states 


probabilities 


comments 


Ql 


a G Sa 


P(a) 




k 


beS b _ 


P(b) 




c 


c G Sc 


P(c\x) 




d 


deSd 


P{d\x) 




X 




P(x\a,b) 





Table 7 
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Figure 11: CB net for 2-body scattering. We show how to construct a q-embeddinj 
for this CB net. 



Applying Step 1 to Af c for two body scattering yields Af c mo d defined by FiglT^l 
and Table |H1 



O b, 



■3 • 




Figure 12: CB net of Fig^2 after adding marginalizer nodes. 
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nodes 


states 


probabilities 


comments 




a 2 e Sg, 


Pa(a 2 ) 






as G Sg, 


<5(a 3 , a 2 ) 




k 2 


b 2 eS b _ 


A(& 2 ) 




h 


b 3 eS b _ 


5(&3,& 2 ) 




0.2 


c 2 e 5c 


-Pc|x(c 2 \x 3c ) 




03 


c 3 e 5"c 


<5(c 3 , c 2 ) 




d 2 


d 2 e Sd 


Pd\x(d 2 \x 3 d) 




d 3 


d 3 e Sd 


5(d 3 ,d 2 ) 




{2L 2 cj 2i2d) 


(x 2 c,x 2 d) e Si 


Px\a,b(x2c\a 3 , b 3 )5(x 2 d, %2c) 




X 3c 


x 3c e 


S(x 3c ,X 2c ) 




X 3 d 


X 3 d G 5s 


S(x 3d ,X2d) 





Table 8 



(Step 2) Replace node probability matrices by their q-embeddings. Add 
ancilla nodes. 

More specifically, replace Af c mo d by a QB net J\f® obtained as follows. For each 
node of Af c m od, except for the marginalizer nodes that were added in the previous 
step, replace its node matrix by a new node matrix which is a q-embedding of the 
original node matrix. Add a new node for each ancilla index of each new node 
matrix. These new nodes will be called ancilla nodes (of either the source or sink 
type) because they correspond to ancilla indices. 

Applying Step 2 to net Af c mo d for two body scattering yields defined by 
FigEHand Table El 
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Figure 13: A QB Net which is a q-embedding for the CB net of FiglTTl 



nodes 


states 


amplitudes 


comments 


SlI 


Ol G So 


<y(oi,o) 




a 2 


»2 G 5a 


A(a 2 \ai = 0) = V p «( fl 2) 






a3 G Sa 


^(03, 02) 






(a,4,b4,X 2c ,X2d) G Sa,b,x,x 


A(o 4 , &4, a;2c, a;2d|a 3 , 63, ari c = 0, arid = 0) = 
\J Px\a,b(x2c\a3,h)S(a i , a 3 )S(b i: b 3 )S(x 2 d, x 2c ) 






a 5 G 5a 


S(a 5 , a 4 ) 




h 


6l G Sb 


tf(6i,0) 






b 2 G 5b 


A(6 2 |&i = 0) = VA(6a) 




h 


63 G S k 


5(63,62) 




h 


65 G 5& 






£1 


Ci G 5 £ 


<5(ci,0) 






(c 2 , £4c) G Sc.x_ 


A(c 2 ,x 4c |ci = 0,x 3c ) = 

^ Pc\x(c2\x 3c )S(x4c, X 3 c) 




£3 


C3 G 5 £ 


S{c 3 ,c 2 ) 






di G 5d 


S(d u 0) 






(^2,^4(i) G Sd,x 


A(d 2 ,ar4d|cii = 0,ar 3rf ) = 
\J p d\x (<k 1 £3d ) 5 {x4d , x 3d ) 




d 3 


d 3 G 


S(d 3 ,d 2 ) 




3hc 


Xi c G 5^ 


S(x lc ,0) 




£ld 


x\d G 5 £ 








X 3c G Sx 


S{x 3c ,X 2c ) 






X3d G Sx 


&{3>3d,X 2 d) 






X§c G 5a; 


5{Xf >c ,X ic ) 






^50" G 5^ 


S(x 5d ,x 4d ) 





Table 9 



looks much more complicated than Af c , but it really isn't, since most of 
its node matrices are delta functions which quickly disappear when adding over node 
states. 

According to Table the probability amplitude for the external (aka leaf) 
nodes is given by 



M a 5, h, c 3 , d 3 , x 5c , x 5d ) = 

= E. \j Pa{ <h) P b{ b 2 ) P x\a,b(x 2c | a 3 , b 3 )Pc\x{c2 I X 3 c) P d\x{d 2 1 X 3d ) 

0(ftj = <h = «4 = a 5)9( b 2 = b 3 = h = h) } ( 115 ) 

9{c 2 = c 3 )9(d 2 = d 3 ) 

9(<h = h = c x = dj_ = x lc = x ld = 0) 

where we have summed over all internal (non-leaf) nodes. Eq. (jll5|) immediately 
reduces to 



A(a 5 b 5 ,c 3 ,d 3 ,x 5c ,x 5d ) = 

= ^/PJosjP^^ = x 5d ) 

Eq. ljllfjjl shows that the net that we constructed from the net M c by following 
steps 1 and 2 satisfies the definition Eq. (jl()5|) that we gave earlier for a q-embedding 
of M c . The probability distribution of the states of the external nodes of the QB net 
contains all the probabilistic information of the original CB net M . Hurray! 

From Eq. (jll6|) . it is clear that by running J\f® on a quantum computer (or 
similar quantum system), we can calculate any conditional probability that one would 
want to calculate for Af c . For example, suppose we wanted to calculate Pa,d\x- R- un 
on the quantum computer several times, each time measuring nodes 05,^3 and 
x 5rf and not measuring all other external nodes. The resulting measurements will 
be distributed according to Pa,d,x- Taking the magnitude squared of the amplitude 
and summing the result over the states of the un-measured external nodes will be 
performed automatically by Nature. The laws of quantum mechanics guarantee it. 
Proceed in the same way to calculate P £ . Run on the quantum computer several 
times, each time measuring node x M and not measuring all other external nodes. 
Finally divide Pg,,d,x by Px on a classical (or quantum?) computer. 

The q-embedding of a CB net, as defined by Eq. (jl(J5|) . is not unique. For 
example, we could have defined the net M Q given by FigHS] without nodes % and b 3 . 
We chose to include such nodes for pedagogical reasons. 

To run a QB net on a quantum computer, we need to replace the QB net 
by an equivalent SEO that a quantum computer can understand. This can be done 
with the help of a quantum compiler |H][E]- One could compile individually each 
node representing a q-embedding, or one could compile whole subgraphs of the QB 
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net all at once. Note that it may suffice to find a SEO that is only approximately 
(within a certain precision) equivalent instead of exactly equivalent to the QB net. 
This may be true if, for example, the probabilities associated with the CB net that 
was q-embedded were not specified too precisely to begin with. 

Suppose ai,a 2 ,...a„ belong to a finite set Sg,, and suppose that they are 
distributed according to a probability distribution P a . What number v of samples aj is 
necessary to estimate Pa within a given precision? This question is directly relevant to 
our method for estimating probabilities by running a QB net on a quantum computer. 
We will not give a detailed answer to this question here. For an answer, the reader 
can consult any book on the mathematical theory of Statistics. An imprecise rule of 
thumb is that if the support of Pa has uq elements, then v must be at least as large 
as z/ ; i.e., one needs at least "one data point per bin" to estimate Pa with any decent 
accuracy. 

We have given a method for calculating, via a quantum computer, the condi- 
tional probabilities associated with a CB net. Does our method have an advantage in 
time complexity with respect to classical methods for calculating the same probabil- 
ities? We will not give a detailed answer to this question here. The answer must be 
yes, sometimes. After all, our method generalizes the algorithms by Deutsch-Jozsa, 
Simon, Grover, etc., and these are known to have a complexity advantage. 

To conclude this section, we will present a second, more complicated example 
of our method of finding a q-embedding for a CB net. A CB net (first given in 
Ref . [T7] ) for lung disease diagnosis is defined by FigH3|and Table ITU1 




Figure 14: CB net for lung disease diagnosis. We show how to construct a q- 
embedding for this CB net. 
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nodes 


states 


probabilities 


comments 


a 


a £ Bool 


P(a = 1) = .01 


Visited Asia? 


b 


b £ Bool 


P(k= l\s= 1) = .60 

7~>/ L 1 1 A\ OA 

P\b = l|s = 0) = .30 


Bronchitis? 


d 


d £ Bool 


P(d = l|e= l,o= lj = .90 

P(d= l|e= 1,6 = 0) = .70 

tj{ j 1 1 . n /-, a\ on 
r\d — 1 e = U, = 1 J = .oU 

11^ n k nA in 
j — 1 e — U, = UJ — .1U 


Dyspnea(trouble breathing)? 


e 


e £ Bool 


P{e\l,t) = %,Z Vi) 


Either TB or Lung Cancer? 


I 


I £ Bool 


E>// 1 1 ^ 1 \ in 
r [L_ = 1|S = 1) = .W 

P(L = i\s = 0) = .01 


Lung Cancer? 


s 


S £ DOOl 


-T (Ji — i-) — -0 


Smokes? 


t 


t £ Bool 


f (* = l|a = 1) = .05 
P(t = l|a = 0) = .01 


Tuberculosis? 


X 


x £ Bool 


^fe= l|e= 1) = .98 
ffe= l]e = 0) = .05 


Positive X-ray? 



Table 10 



If one follows the two steps that were described earlier in this section, one 
obtains the QB net defined by FigH^l and Table [TT] 
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nodes 


states 


amplitudes 


— i 


ai £ Boo/ 


5(ai,0) 

^ j_ 7 w / 


—2 


«2 G Bool 


A(a 2 K = 0) = A /B a (a2) 




a 3 G Boo/ 


8(a 3 ,a 2 ) 




(a 4 ,i 2 ) G Bool 2 


A(a4,t 2 \a 3 ,t 1 = 0) = i/Ptiaih a 3 )5(a 4 ,a 3 ) 
V -'- 




&5 G Boo? 


<5(a 5 ,a 4 ) 


kl 


bi G Boo/ 


^i,0) 




(&2, S4b) G Bool 2 


A(b 2 ,S4b\bl =0, S 3fc ) = ^/Pt| £ (&2|s3b)'5(s4b,S3b) 


ft, 

—O 


b 3 G Boo/ 


^(63,^2) 


m. (7„ p, ,] 
Vil4 ' —2 5 ^Ad) 


(bA d't eArt) G Bool 3 


A(bA dn f>Art xb-x d,i — pq^l = * / (d,o\bri pr>,i)()(b\ biifilfAj p^^] 




O5 t: DOOL 


0(0 5 , O4) 


Ml 


1 7 -1 T^rtrt] 
(-1 1 t JDUUL 




—3 


d 3 G Bool 


<5 fcf s . d?) 


61 J 

—la 


e\d G Bool 






ei x G Bool 


5(e-\r 0) 




{&2di e 2x, h, ti) 


A( e 2d, e 2 x, ^4, *4 \eu = 0, e la; = 0, Z 3 , t 3 ) = 




G Boo/ 4 


\J P e\Lt{e2d\h,t 3 )S(e2x, e 2 d)6{hi l 3 )8(U, h) 


—3d 


e 3rf G Boo/ 


S(e 3 d, e?d) 


■^3x 


e3a; G Boo/ 


5(e 3x , e2x) 


(jLAx 1 £.2 ) 


(e4x, 2:2) G Boo/ 2 


A(e ix ,x 2 \e 3x ,xi = 0) = y / B ; r|e(^2|e 3a; )(5(e4 a; , e 3x ) 


— bd 


esrf G Bool 


<5 f e^r/, e4r/) 


—bx 


&5x G Boo/ 


<5(e^ T . e 4 r) 


—i 


/1 G Boo/ 




(Zn S*l) 


(7, S4/ 1 ) G Bool 2 


Ado sm\1-\ =0 su) = * Pu , (li\sti)5( sai sv) 


7 

i3 


63 t SDUOi 


0(1-3, h) 


7 

i5 


£5 G -DOOL 


Atl 1 \ 


£l& 


s\i) G Bool 


o(sib, V) 




S\i G Bool 


o(su,v) 




^26, 52/ J G -DOOt 


^4(s2b, S2i|sib = 0, Sll = 0) = -y' Ps_{s2b)5{s2l , S2b) 


^3b 


^36 ^ Bool 


0~(s 3 b, &2b) 




S31 G ±3001 


0(S 3 l, S2l) 


i!5b 


555 G Bool 


0~(s5b, S4b) 




55/ G Bool 


S(S5i, S41) 




fi C Jinn] 


0\li, V) 


*3 


t 3 G Boo/ 


S(t 3 ,t 2 ) 


& 


£5 G Boo/ 


6{t B ,U) 




Xi G Boo/ 


S(x u 0) 


£3 


x 3 G Boo/ 


S(x 3 ,x 2 ) 



Table 11 
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According to Table \H\ the probability amplitude for the external (aka leaf) 
nodes is given by 



A( a 5, h, d 3 , e 5d , e 5x , l 5 , s 5b , s 5l , t 5 , x 3 ) 



Pg£a 5 )Pt\a(t 5 \a 5 )Pb\s(b 5 \s 5b )P^^^ 



0(e 5d ,e 5x )6(s 5b ,s 5 i) 



(117) 



6 Voting Net and Grover's Microscope 

In this section we will first present a CB net, call it M c ', that describes voting. 
Then we will find a QB net that is a q-embedding of M c . In certain cases, the 
probabilities that we wish to find are too small to be measurable by running on 
a quantum computer. However, we will show that sometimes it is possible to define a 
new QB net, call it J\f® , that magnifies and makes measurable the probabilities that 
were unmeasurable using alone. We will refer to as Grover's Microscope 
for J\f® , because J\[® is closely related to Grover's algorithm, and it magnifies the 
probabilities found with N®. 

Suppose y G Bool and x = (x ,^ 1 , . . . ,1^"') G BooI Nb . Consider the CB 
net ( "voting net" ) defined by FigH^] and Table IT21 




Figure 16: "Voting" CB net. 



nodes 


states 


probabilities 


comments 


for all % G Z 0i at b -i 


x l G Bool 


P{x i ) 




y 


y G Bool 


P(y\x) 





Table 12 



Henceforth, we will abbreviate P(y = 0\x) = Pi and P(y = l\x) = qi, where 
i = dec(x) G Zo,jv s -i- Hence Pi + qi = 1 for all i G Z 0jNs _ 1 . In general, the probability 
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matrix P(y\x) has 2 Nb free parameters (namely, Pi for all i G Zq^ s _ 1 ). This number 
of parameters is a forbiddingly large for large Nb- To ease the task of specifying 
P(y\x) , it is common to impose additional constraints on P(y\x). An interesting 
special type of P(y\x) is deterministic pd(Bool\Bool NB ) matrices; that is, those 
that can be expressed in the form 



P(y\x)=5(yJ(x)), (118) 

where / : BooI Nb — > Bool. In this case, the voting net can be used to pose the 
satisfiability problem (SAT): given y = 0, find the most likely x G BooI Nb ; in 
other words, find those x for which f(x) = 0. We say / is AND-like if all Pi equal 
zero except for one pi which equals one. For example, for Nb — 2, if / is an AND 
gate, then 



(x°, x l ) 



P(y\x] 



AND 





00 


01 


10 


11 


yi o 


1 


1 


1 





l 











1 



(119) 



matrices; that is, those that can be expressed in the form 

p(vW) = E *(v, f^))P{t>")P{t l \x l ) ...Pit 



-U X N B - 



(120) 



where / : BooI Nb Bool and we sum over all t = (t°, t 1 , . . . , t^" 1 ) G BooI Nb . 
When f(t) = Pvt 1 V. . .Vt NB ~ l , P(y\x) is called a noisy-OR. Appendix lAl discusses 
how to q-embed deterministic pd(Bool\Bool NB ) matrices, and how to express such q- 
embeddings as a SEO . Appendix El discusses the same thing for quasi-deterministic 
pd{Bool\Bool N}3 ) matrices. 

A q-embedding for the CB net defined by Fig^Jand Table IT21 is given by the 
QB net defined by FigEI and Table ED 



nodes 


states 


amplitudes 


comments 




xi G BooI Nb 


5(xx,0) 






x 2 g BooI Nb 


A(x 2 \x! =0) = ^/Ps(x 2 ) 






(£3,2/2) G Bool N * +1 


A(x 3 ,y 2 \x 2 ,y! = 0) = ^IPy^(y 2 \x 2 )S(x3,x 2 ) 




zu 


x 4 G BooI Nb 


5(X4,X3) 




v, 


yi G Bool 


S(yi,o) 






y 3 G Bool 


<%3,J/2) 





Table 13 



According to Table EH the probability amplitude for the leaf (external) nodes 



is 
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Figure 17: A QB Net which is a q-embedding for the CB net of FigfTBI 



Yl ^ p i^2) p y\x{y 2 \x 2 )6(y 2 = y 3 )6(x 2 = x 3 = £ 4 )#(?i = y 1 = 0) (121a) 




To fully specify the QB net for voting, we need to extend A(x*2\xi = 0) and 
^4(^3; 2/2^2) yi — 0) into unitary matrices by adding columns to them. This can always 
be accomplished by applying the Gram-Schmidt algorithm. But sometimes one can 
guess a matrix extension and applying Gram-Schmidt becomes unnecessary. If Pg is 
uniform (i.e., P(x) = 1/Ng for all x, which means there is no a priori information 
about x), then A(x 2 |£i = 0) = 1/y/Ns- In this case, we can extend A(x 2 \xi = 0) into 
the unitary matrix 



(This works because all entries of the first column of Hn b are equal to l/y/Ng-) As 

to extending A(xs, y 2 |^2, yi = 0), this can be done as follows. Define 



[y4(f 2 |£l)] = H Nb . 



(122) 





(123) 



and 




(124) 



A possible way of extending A(xs, 2/2^2, 2/i = 0) into a unitary matrix is 




(125) 
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Unitary matrices of this kind are called D- matrices in Ref. [H]. Ref. [Sj shows how to 
decompose any D-matrix into a SEO. 

Earlier, we explained how to estimate a conditional probability for a CB net by 
running a QB net v times on a quantum computer. If we wanted to find P{y\x°,x 1 ) 
for the voting CB net, then the number of runs v required to estimate Piy^^x 1 ) 
with moderate accuracy would not be too onerous, because the domain of P(y\x°, x 1 ) 
is Bool 3 , which contains only 8 points. But what if we wanted to estimate P(y\x)7 
For large Nb, the domain of P(y\x) is very large (2 Nb+1 points). If the support of 
P(y\x) occupies a large fraction of this domain, then the number of runs v required 
to estimate P{y\x) with moderate accuracy is forbiddingly large. However, there are 
some cases in which "Grover's Microscope" can come to the rescue, by allowing us to 
amplify certain salient features of P{y\x) so that they become measurable in only a 
few runs. 

Next we will discuss Grover's Microscope for the voting QB net defined by 
FigHZI and Table For simplicity, we will assume that P% is uniform. 

Let k = (k , Kjv s -i) label Nb bits and let r label another bit. Assume 

that t and all the Kj £1X6 distinct. Define 



)k = 4> P = (y/po, VPii • ■ ■ > y/PN s -i) 



(126) 



and 



HI K 



Qo, VQi 



(127) 



I*) 



1 



s 



(|0) r |^) s +|l) T |^)«) 






:i28i 



Since pi + qi = 1 for all i, (j)p(j)p + (j) q (j)q = N$- According to Eq. (jl21|) . when P$ is 
uniform, the voting QB net fully specifies a unitary matrix U net such that 



|¥) = EU|0)*|0), 
Define orthonormal vectors eo and e± by 



(129) 



e 







ei 



(130) 



where V is a unit vector in the direction of V. If P(y\x) is deterministic with AND-like 
/, then all components of eo are zero except for the one at the target state jtarg- 
In terms of eo, ei, ^ can be expressed as 



11 



1 



■■{\4> P \eo + \4> q \ei 



It is convenient to define a vector orthogonal to 

1 



: (l0 Q |eo - |0 P |ei) • 



(131) 



(132) 



If P(|/|x) is deterministic with AND-like /, then 



1 and 



large N s , \& ~ ei and 

*<,= 1 



eo- For an arbitrary angle a, let 
c«|0,| +s-|0 p |)e o + -c°|0 p |)ei 



V^s — 1 so, for 



(133) 



where sa = sin A and ca = cos A for any angle A. Let l(x,y) denote the angle 
between 2 vectors x and y. Note that Z(*j_, = a/2. We define Z(ei, = 0/2. 

Fig^J portrays various vectors that arise in explaining Grover's Microscope. 
Note that = e$ when a = 0. 




Figure 18: Various vectors relevant to Grover's Microscope. 



Since we plan to stay within the two dimensional vector space with orthonor- 
mal basis eo, ei, it is convenient to switch matrix representations. Within span(eo, ei), 
eo, ei can be represented more simply by: 

e ° = ( o ) ' 61 = ( i ) • (134) 

If eo, ei are represented in this way, then 




42 



and 




■qf' ± = WV , where W = 7 'in- ( 137 ) 



The matrix y ^ J J is a clockwise rotation by 7r/2 in space span(eo, e\). Thus, 

equals a clockwise rotation by tt/2 followed by a counter-clockwise rotation by a/2. 
Define the following reflection operators 

Ro = l- 2U ]qU U 1q)t = (-l)%>* n i<»r } (138) 

RiS, = U ne tRoUl et , (139) 

Rv ± = WR*W* = WU net R Ul t W^ . (140) 
From Eq.(|2I|), it follows that 

- RyRy^ = c a ^ T - s a Wl + s a ^ ± ty T + c a y ± ^l . (141) 

Thus, —R^R^,' ± rotates vectors in span(eo,ei), clockwise by an angle a. 

Grover's Microscope can be summarized by the following equation 

(-i^/J'-tf « eo , (142) 

for some integer r to be determined, where "~" means approximation at large N$. 
What this means is that our system starts in state \& and is rotated consecutively 
r times, each time by a small angle a, until it arrives at the state e$. If P(y\x) is 
deterministic with AND-like /, then measuring state eo yields the target state jtarg- 
The optimum number r of iterations is 

ra«|(l + 2fc) (143) 

for some integer k. Note that cos(6 l /2) = (W|ei) = \4> g \/\/Ns so, in general, 9 depends 
on \(j) p \ (or on \<f> q \ = J 'Ng — \4> p \ 2 )- If P(y\x) is deterministic with AND-like /, then 
\<pp\ = 1 and \(j) q \ — y/N$ — 1. In this case, it is convenient to choose a = 8, so that 
ty' ± = e\ and FigsJU] and [IE] become the same diagram under the mapping ^ — > \l 
and — > <p = Cq. Then the optimum number r of iterations for Grover's original 
algorithm and for Grover's Microscope are equal. If we don't know ahead of time the 
value of \4> p \, then setting 9 = a will make both r and a depend on the unknown \(f) p \, 
although the product ra will still be independent of it. 
Let 
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V fiscope 




-e±e% + e ef 



= -W T L + . (144) 

Note that 



U^ scope ^ = *1 . (145) 

From the point of view of quantum compiling, Grover's Microscope approximates 
the 7r/2 rotation U^ scope by the r-fold product of —Rq,R^>^ where we assume that 
—Rq,Rq,' ± can be shown to have a SEO of low (polynomial in Nb) complexity. (If 
such a low complexity SEO cannot be found, then it is pointless to divide U^ scope into 
r iterations of —R^R^,' , and we might be better off compiling U^ scope all at once.) 



A Appendix: Deterministic pd(Bool\ Bool B ) matri- 
ces 

In this Appendix, we will first define a special kind of probability matrices which we 
call deterministic pd(Bool\Bool NB ) matrices. Then we will show how such probability 
matrices can be q-embedded, and how their q-embedding can be expressed as a SEO. 
Suppose y G Bool and x — (x ■ X j • • • 2 x Nb ^) G BooI Nb . Let / : BooI Nb -> 

Bool. 

We will say that / is AND-like if f(x) = 8(x = x tar g) for some target vector 
%targ £ BooI Nb . An AND-like / maps all x into zero except for x tar g which it maps 
into one. Thus, |/ _1 (1)| = 1. An example of an AND-like / is the multiple AND gate 
f(x) = x° Ax 1 A. . .Ai*" 1 , which can also be expressed as f(x) = 9[x = (1, 1, ... , 1)]. 

We will say that / is OR-like if f(x) = 9{x ^ Xtarg) for some target vector 
Xtarg £ BooI Nb . An OR-like / maps all x into one except for Xtarg which it maps 
into zero. Thus, |/ _1 (0)| = 1. An example of an OR-like / is the multiple OR gate 
f(x) = x°Vx 1 V. . .Vi^ 8 " 1 , which can also be expressed as f(x) = 8[x ^ (0, 0, ... , 0)]. 

We will say that / has a single target if it is either AND-like or OR-like. If 
/ has more than one target (i.e., if |/ _1 (0)| and |/ _1 (1)| are both greater than one), 
then we will say that / has multiple targets. 

Suppose y G Bool and x = (x ,^ 1 , . . . .x 1 ^ 13 ^ 1 ) G BooI Nb . Let / : BooI Nb — * 
Bool. In this section, we consider deterministic pd(Bool\Bool NB ) matrices; that 
is, probability matrices of the form P(y\x) = 5(y, f(x)). First let us consider the case 
that / has a single target. For example, for Nb = 2, if / is an AND gate 
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P(y\x] 



AND 



and if / is an OR gate 





00 


01 


10 


11 


yi o 
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all j e Z { 



Suppose bit value y is stored in the bit labelled r 



x 



0,N S - 



And suppose bit values 
are stored in the bits labelled k = (k , Ki, . . . , kn b -i)- Define e,- for 
i to be the N$ dimensional column vector with jth component equal 



ejej and Il targ 



n 



Jtarg 5 



to one and all other components equal to zero. Let Uj 
where j tar g £ ^o,jv s -i is the target state. II tar5 can expressed as product of number 
operators. Indeed, if 



Jtarg 
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^ ] Xtarg,i^ 
i=0 
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then 



N B -1 



^-targ ~ ^-jtarg 



[n(Ki)6(x t arg,i = 1) + n(Ki)9(x ta rg,i = 0)] 



(149) 



i=0 



For example, if j targ = then H targ = n(K M«i) . . . ra(«jv B _i). 

An AND-like probability matrix P(y\x) is q-embedded within the unitary ma- 
trix 



U 



AND- 





y = 


y = l 


y = 





1 — n ta rg 


-fftarg 


y = 


1 


^-targ 
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Note that 



(150) 



U AND- 



like 



-1 -1 
1 -1 



n 



targ 



1 + U targ (K)(-ia y (T) - 1) 
[-zt7 y (r)] ntar «^ . 



(151a) 

(151b) 
(151c) 



Eqs. (ll49|) and ()151cj) show how to express U^ND-Uke as a qubit rotation with multiple 
control qubits. Operations of this kind can be decomposed into a SEO using the 
techniques of Refs. jlj and [Bj. 
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An OR-like probability matrix P(y\x) is q-embedded within the unitary matrix 

(152) 
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Note that 
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(153a) 
(153b) 



Finally, let us consider the case when / : BooI Nb — > Bool has multiple targets. 
Let T C Z QtNs _i be the set of these targets; i.e., either T = / _1 (0) or T = / _1 (1). 
Define n targ by 



n 



targ 



J6T 



(154) 



ritarg can be expressed as a product of number operators. Indeed, each Ylj on the 
right hand side of Eq. (jl54j) can be separately expressed, using Eq. (|149|) . as a product 
of number operators. If T = /~ 1 (1), then P(y\x) is q-embedded within the unitary 
matrix 



multi—targ 



[A(y,x\y,x)] 
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(155) 



B Appendix: Quasi-deterministic pd(Bool\Bool B ) 
matrices 

In this Appendix, we will first define a special kind of probability matrices which 
we call quasi-deterministic pd(Bool\Bool NB ) matrices. Then we will show how such 
probability matrices can be q-embedded, and how their q-embedding can be expressed 
as a SEO. 

Suppose y G Bool and x = (x ,^ 1 , . . . ,x Ng ^ 1 ) G BooI Nb . Let / : BooI Nb — > 
Bool. In the previous appendix, we considered deterministic pd(Bool\Bool NB ) matri- 
ces; that is, probability matrices of the form P(y\x) = 5(y, f{x)). In this section, we 
will consider quasi-deterministic pd(Bool\Bool NB ) matrices; that is, probability 
matrices of the form 
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P(y\x) =J2s(yJ(t))P(t°\x°)P(t 1 \x 1 ) . . .p^- 1 !^- 1 ) 



(156) 



where we sum over alH = (t°, t 1 , . . . , t NB *) £ BooI Nb . FiglT^l shows a CB net repre- 
sentation of Eq. (J156|) . Examples of quasi-deterministic pd(Bool\Bool NB ) matrices are: 



X 



X 



X 



X N B- 1 




Figure 19: Quasi-deterministic ("noisy") pd(Bool\Bool NB ) gate. 



(l)the noisy OR, for which f\t) = t° Vt 1 V . . . Vt^ 5 " 1 ; (2)the noisy AND, for which 
f(t) =t°At 1 A...At NB - 1 ; (3)the noisy CNOT, for which fit) = t° m 1 ® . . .®t NB ~\ 
etc. 

For each a £ Z 0tNg _i, the probabilities P(t a = = x) will be abbreviated 
by pf x for t,x £ 5ooZ. P(t a = = x) has two independent parameters which we 
may take to be Pqi (the probability of false negatives) and Pi Q (the probability of false 
positives). Pq and p™ x can be expressed in terms of these independent parameters: 
Poo = 1 — Pi , Piy = 1 — Poi- Whereas a completely general probability matrix 
P{y\x) £ pd{Bool\Bool NB ) has 2 Nb free parameters, a quasi-deterministic P(y\x) has 
2Nb free parameters. 

Rather than q-embedding the probability matrix P(y\x) as a whole, it is con- 
venient to q-embed separately the probability matrices P(y\t) and P(t a \x a ) for every 
a £ Z 0tNg _i. P(y\t) = 5(y,f(i)) is a deterministic pd(Bool\Bool NB ) matrix so its q- 
embedding is discussed in Appendix |XJ As for P(t a \x a ), it can be easily q-embedded 
as follows. For each a £ Z ^ Nb _i, let 

Jk)' sk)- (157) 

P(t a \x a ) is q-embedded within the unitary matrix: 

[A(t a ,S a \t a ,x a )}= ( ^ ) . (158) 

Unitary matrices of this kind are called D- matrices in Ref. [Hj. Ref. [B] shows how to 
decompose any D-matrix into a SEO. 
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